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We theoretically study the crossover between spin Hall effect and spin swapping, a recently pre¬ 
dicted phenomenon that consists in the interchange between the current flow and its spin polarization 
directions [Lifshits and D’yakonov, Phys. Rev. Lett. 103 , 186601 (2009)]. Using a tight-binding 
model with spin-orbit coupled disorder, spin Hall effect, spin relaxation and spin swapping are 
treated on equal footing. We demonstrate that spin Hall effect and spin swapping present very 
different dependences as a function of the spin-orbit coupling and disorder strengths. As a conse¬ 
quence, we show that spin swapping may even exceed spin Hall effect. Three set-ups are proposed 
for the experimental observation of the spin swapping effect in metals. 


I. INTRODUCTION 

The spin-orbit coupling (SOC) is a relativistic effect 
that couples the particle’s spin degree of freedom to its 
orbital angular momentum. In solid-state, this coupling 
originates from the interaction between the carrier’s spin 
and the magnetic field that it experiences in its rest 
frame in the presence of a potential gradient (crystal 
field, defects etc.). The locking between spin and or¬ 
bital angular momenta has spectacular consequences in 
normal metals such as, but not limited to, spin Hall 
effect spin galvanic effectand spin relaxation^^^. 
In semiconductors^^, metals^^”^^ and more recently topo¬ 
logical insulators^^, SOC is now commonly engineered 
and exploited to generate pure spin currents, thereby en¬ 
abling the electrical manipulation of the spin degree of 
freedom in the absence of an external magnetic field. 

The most emblematic effect induced by SOC is prob¬ 
ably the spin Hall effect (SHE) originally predicted by 
Dyakonov and Perel^ and revived thirty years later by 
Hirsch^^ and Zhang^^. In analogy with the ordinary Hall 
effect, SHE consists in the generation of a transverse spin 
current in response to an electric field applied in the lon¬ 
gitudinal direction. Geometrically, the spin current has 
the form Jj = e^jkOLn^cjo/e)Ek^ where Jj is a spin cur¬ 
rent polarized along the direction i and flowing along the 
direction j, Ek is the electric field applied along the direc¬ 
tion k, (Jo is the longitudinal conductivity, q^h is the spin 
Hall angle (measured in %) and Cijk is Kronecker symbol. 
In their pioneering theory, Dyakonov and Perel derived 
SHE from Mott scattering on spin-orbit coupled impuri¬ 
ties, i.e. from the asymmetric spin-dependent scattering 
of an initially unpolarized electron flow^^’^^. This mech¬ 
anism is usually referred to as extrinsic since it depends 
upon the presence of impurities in the system. Another 
extrinsic effect that occurs upon scattering is the shift 
of the position of the incoming wave packet, called side 
jump scattering^^. This shift renormalizes the velocity 
operator by creating an effective anomalous velocity pro¬ 
portional to the number of impurities in the system^^. 


Besides extrinsic contributions to SHE, the spin-orbit 
coupled band structure may also dramatically influence 
the velocity of traveling electron wave packets even in the 
absence of impurity scattering. Indeed the velocity oper¬ 
ator, V = dpH, is affected by the SOC present in the band 
structure which induces an additional Lorentz force that 
bends the trajectory of Bloch states^^. The effective mag¬ 
netic field felt by the moving electron (also called Berry 
curvature^^’^^) and arising from the band structure is in¬ 
trinsic as it does not rely on impurity scattering. This 
contribution has been widely studied in the context of 
anomalous Hall effect^^, for instance. The first prediction 
of impurity-free intrinsic SHE in non-centrosymmetry 
semiconductors^^’^^ paved the way towards the predic¬ 
tion and discovery of topological insulators^^. 

The search for efficient SHE has been quite intense in 
the past ten years^^~^^. After its pioneering observation 
in semiconductors^’^, the attention has quickly drifted 
towards the exploration of metals^^. A wide variety of 
noble and transition metals has been investigated and 
while the exact values and the proper method to detect 
SHE is still under debate, clear indication of intrinsic 
SHE^^’^^ has been established in Ad and bd noble metals 
(Pt, Pd, Ta, W etc. - see Ref. 30-35), as well as in M 
transition metals^^. Alternatively, it has been recently 
suggested that large spin Hall angle q^h could be achieved 
by exploiting skew scattering on the resonant states of 
bd or 4/ impurities embedded in a light metal matrix, 
such as Cu or Al^^”^^. Evidence of large Hall angles 
from resonant skew scattering was reported in Pt-^^’^^ 
and W-doped Au^^, as well as Ir-^^ and Bi-doped Cu^^. 

Metals doped with heavy impurities present an inter¬ 
esting paradigm as they may be utilized to observe novel 
effects such as spin swapping (SSW), i.e. the conversion 
of a spin current Jj into a spin current J/ upon scatter¬ 
ing by spin-orbit coupled impurities^^’^^. In the present 
work, we explore the nature of spin-orbit coupled trans¬ 
port within a disordered tight-binding model and demon¬ 
strate the emergence of SSW when the injected current 
is spin-polarized. In our model, SSW, SHE and spin 
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relaxation all arise from the same spin-orbit coupled dis¬ 
order which allows us to treat these different effects on 
equal footing. In particular, a crossover between SSW 
and SHE, controlled by the spin relaxation, is revealed 
and analyzed. This article is organized as follows. The 
physics of impurity-induced spin swapping is first out¬ 
lined in Section II. In Section III, a single band tight- 
binding model and the corresponding non-equilibrium 
Green’s function technique are presented. Section IV 
is devoted to the numerical results. Finally, Section V 
suggests three setups for the experimental observation of 
SSW, and Section VI concludes this study. 


II. PHYSICS OF SPIN SWAPPING 
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FIG. 1. (Color online) (a) Sketch of Mott scattering by a spin- 
orbit coupled impurity. Electrons with an out-of-plane spin 
polarization pointing up have a larger probability to scatter 
towards the right while electrons with an out-of-plane spin 
polarization pointing down have a larger probability to scatter 
towards the left, (b) Sketch of spin swapping mediated by a 
spin-orbit coupled impurity. When the spin of the carrier lies 
in the scattering plane (fc, k'), it experiences a magnetic field 
Bso oc k' xk that depends on the direction towards which the 
spin is scattered. This induces a scattering-induced precession 
giving rise to the spin swapping. 

The physics of extrinsic Mott skew scattering and SSW 
is illustrated on Fig. 1(a) and (b), respectively. Let us 
first consider spin-dependent scattering against a spin- 
orbit coupled impurity. In real space, the impurity po¬ 
tential reads 

^ Vimp{r - n) + {Xso/h)&■ (VMmp(r - n) X p), 

i 

( 1 ) 

where Vmp is the spin-independent impurity potential, 
and Ago is the SOG parameter of the impurity. In the 
reciprocal space, this impurity potential becomes 

HZ' = (2) 

i 

It is clear that the spin-orbit coupled part of the impu¬ 
rity potential acts like a magnetic field oc k' x k 
on the incoming electron spin (T, where k is the momen¬ 
tum of the incoming electron and k' is the momentum 


of the outgoing electron. Therefore, in the case of an 
unpolarized charge current, this magnetic field defines a 
local quantization axis so that one can phenomenologi¬ 
cally separate electrons with a spin polarization pointing 
parallel and antiparallel to Hgo, as illustrated on Fig. 
1(a). In this case, proper treatment of the collision inte¬ 
gral leads to Mott scattering, i.e. electrons with a spin 
momentum pointing (anti)parallel to k' xk have the ten¬ 
dency to scatter towards the left (right). The resulting 
spin current is therefore polarized along k' xk and flows 
transversely to both k' xk and k. This mechanism is at 
the core of the skew scattering mechanism giving rise to 
SHE as originally pointed out by Dyakonov and Pered. 

Let us now consider an incoming spin-polarized current 
(it can be either a pure spin current or accompanied by a 
charge flow), whose spin polarization lies in the scattering 
plane {k^k'). In this case, the incoming spin precesses 
around the effective magnetic field B^q. Since this mag¬ 
netic field only exists upon scattering from k to k\ the 
spin polarization of the outgoing current is re-oriented^^. 
Lifshits and Dyakonov tagged this effect ” spin swapping” 
as the spin polarization and flow direction of the incom¬ 
ing spin current are swapped during this process: an in¬ 
coming spin current Jj gives rise to a spin current J- 
when i 7^ j (similarly a spin current, say, produces 
two spin currents and J|). 

Physically, SSW is nothing but a spin precession 
around the spin-orbit field Hgo- Such precessions are well 
known in semiconductors possessing Rashba, Dresselhaus 
or Kohn-Luttinger SOG^^. This SOC-induced preces¬ 
sion has been utilized in the Datta-Das transistor^^’^^ 
(see also Ref. 47) and its impact on spin transfer torque 
has been investigated in (Ga,Mn)As-based spin-valves^^ 
and single layers^^. However, the impurity-induced spin 
swapping presents a number of interesting differences. 
First, in contrast with coherent spin precession around 
Rashba or Dresselhaus SOG, spin swapping arises from 
spin precession stemming from incoherent scattering. It 
is quite remarkable that this effect survives the disorder 
configurational average, as does extrinsic SHE. A spin 
diffusion equation in the presence of spin-orbit coupled 
impurities within the first Born approximation has been 
derived by Shchelushkin and Brataas^^ a few years be¬ 
fore the prediction of SSW by Lifshits and Dyakonov. 
The charge and spin current equations obtained in Ref. 
52 are reproduced below 

jc/cro = -V/ic + X ps, (3) 

eJi/ao = -(1 + X Vmc + 

(4) 

Here, jj^c is the spin-independent chemical potential, ii\ 
is the i-th component of the spin-dependent chemical po¬ 
tential, and ^so = Aso^f unitless spin-orbit param¬ 

eter. The metallic system is described in terms of its free 
electron Fermi wave vector and its mean free path A. 
jc is the charge current density and Ji is the spin current 
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density defined as a vector in the direction of the spin po¬ 
larization, flowing along the direction e^. The first terms 
in Eqs. (3) and (4) are the diffusion terms, the second 
terms (oc are the side jump contribution producing 
SHE [Eq. (4)] and inverse SHE [Eq. (3)]. The third term 
in Eq. (4) is the spin swapping effect. Note that since 
these equations are derived within the first Born approx¬ 
imation (i.e. up to the second order in impurity potential 
only), skew scattering is neglected (see also Ref. 
53). The spin swapping term in Eq. (4) clearly con¬ 
verts into Ji, which is consistent with Lifshits and 

Dyakonov theory^^. A simplified version of the drift- 
diffusion equations, Eqs. (3)-(4), has been numerically 
investigated by Sadjina et al.^^, neglecting SHE. 

A second interesting aspect of SSW is that since it 
comes from impurity scattering, it can be engineered 
by using resonant scattering on heavy impurities em¬ 
bedded in a light metal host as suggested previously for 
SHE^^”^^. This aspect is beyond the scope of the present 
work, but deserves further investigations. 


III. TIGHT BINDING MODEL 


spin accumulation is relaxed, over a distance of the order 
of the spin diffusion length Agf, SHE smoothly takes over 
and a spin accumulation at the edges of the sample and 
polarized along the normal of the plane emerges. We 
numerically demonstrate this scenario using the tight- 
binding model described below. 
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The drift-diffusion model presented in the previous sec¬ 
tion, Eqs. (3)-(4), has been obtained within the first Born 
approximation, i.e. only accounting for effects propor¬ 
tional to the impurity density and up to the second order 
in impurity potential^^’^^. Therefore, configurational av¬ 
eraging is performed and skew scattering is neglected. In 
this model, the ratio between SHE [second term in Eq. 
(4)] and (side-jump only) SSW [third term in Eq. (4)] is 
oc 1/ AA^f- 

The objective of the present work is to investigate 
the nature of spin-orbit coupled transport (i.e. SSW, 
SHE and spin relaxation) over a broad range of disorder 
strengths and accounting for all the relevant contribu¬ 
tions to extrinsic SHE (side-jump, skew scattering etc.). 
To do so, we consider a tight-binding model on which 
spin-orbit coupled disorder is implemented. By comput¬ 
ing the non-equilibrium local spin density that accumu¬ 
lates at the edges of the sample, one can evaluate the 
relative magnitude of SSW and SHE, identify its depen¬ 
dence as a function of the disorder strength and clarify 
the role of spin relaxation. 

The metallic system we consider is sketched in Eig. 
2(a) and (b). It is composed of a ferromagnet (EM) po¬ 
larized along y and a normal metal (NM) with spin-orbit 
coupled disorder (see below). Eollowing the discussion 
given in the previous section, we choose the EM polar¬ 
ization along y to insure that the spin density injected 
in the normal metal does not readily undergo SHE as 
long as the spin polarization is conserved [see Eig. 2(a)]. 
Therefore, we expect SSW to take place in a region close 
to the interface and limited by the spin relaxation length. 
This SSW generates a spin accumulation at the edges of 
the sample whose polarization is aligned with the direc¬ 
tion of injection (x in this example). While the injected 


FIG. 2. (Color online) (a) An illustration of the crossover 
between SSW and SHE: SSW produces an edge spin accumu¬ 
lation that is polarized along the direction of injection (green 
arrows) and only survives close to the interface, while SHE 
builds up smoothly and results in a spin accumulation with a 
polarization normal to the plane, (b) Schematic of the tight 
binding model of the system. The circles refer to the atomic 
sites and e stands for the on-site energy. The right and left 
leads are free from disorder with e = eo- The central region 
is composed of a ferromagnet with polarization parallel to y 
and the disordered non-magnetic layer (symbolized by atoms 
having different on-site energies eij = eo + , ^ij being the 

random onsite potential). The disorder-induced SOC is sym¬ 
bolized by the fc-dependent magnetic field B{k^k'). The red 
arrows in the non-magnetic region stand for the carriers spins 
which we chose to have random directions to emphasize the 
effect of their coupling with the deflected momenta. 

The calculations carried out in the present work are 
based on the non-equilibrium Green’s function technique 
implemented on a single band tight-binding model^^. All 
the expectation values of the physical quantities of inter¬ 
est are given in real space. The whole system sketched in 
Eig. 2(b), i.e. the EM/NM bilayer, is connected to two 
non-magnetic semi-infinite leads whose role boils down 
to maintaining the flow of the particles throughout the 
system and hence promoting the non-equilibrium regime. 
The total Hamiltonian of the central system reads 

Hq ^ ^ {(^ijf^crcr' T h.C.^ 

~ + + k.C.) -f- Hgo^ 


( 5 ) 
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Here the first term at the right-hand side of Eq. (5) is the 
spin-independent onsite energy in which eij = eo eo 

being the onsite energy constant and being a random 
onsite energy that introduces disorder in the system. The 
local disorder strength is randomly chosen such that G 
[—r/2,r/2]. The second term stands for the exchange 
interaction between the spin of the carriers and the local 
magnetic momentum of site (i, j), is the exchange 
splitting, is the unit vector along the direction of 
the local magnetic moment on site (i, j). In the present 
case, = y in the FM layer and Aij = 0 in the NM 
region. The operator clj ^ creates (annihilates) 

a particle with spin a at position (i, j). The third term 
in the Hamiltonian corresponds to the nearest neighbor 
hopping energy parameterized by the hopping integral t. 

The SOC Hamiltonian Hgo is given by the second 
term in Eq. (1) and can be discretized on the square 
lattice we consider. Since two spatial gradients are in¬ 
volved (VEimp and p = — ifiV), this term yields a spin- 
dependent next-nearest neighbor hopping contribution 
and can be rewritten^^ 


Hso = -i{^so/a^)e{ra){ra x n) ■ (6) 


where are unit vectors pointing from the initial posi¬ 
tion to the nearest neighbor and from the nearest neigh¬ 
bor to the next-nearest neighbor, respectively, and e{ra) 
is the onsite energy of the nearest neighbor. Hence for 
a two-dimensional system, Hgo reduces to two kinds of 
hopping involving the disordered-onsite energies^^ 

where a = \so/o? and a is the lattice parameter. The 
hopping from the site (z — 1, j — 1) to the site {i^j) can 
either pass by (i, j — 1) or by (i — 1, j) which give opposite 
contributions to the spin-orbit field. The same reasoning 
applies for the hopping from site (i — 1, j) to site (i, j — 1). 
In the present study, we choose the tight-binding param¬ 
eters to be in units of the hopping energy t: the transport 
energy E = the exchange interaction A between the 
carrier spin and the local magnetic moment is, unless 
stated otherwise, given by A = t. The SOC and disorder 
strengths {a and T) range from (a = 0.2to(a = 0.6 and 
T = 0.2 X t to T = t respectively). 

The effects we are interested in - SSW, SHE and spin 
relaxation - arise from the conjunction of disorder and 
SOC and hence only emerge when configurational aver¬ 
aging is properly carried out. The advantage of the tight- 
binding approach is that once the calculations are con¬ 
verged, mechanisms such as skew-scattering^^’and side 
jump^^ are fully accounted for. Therefore in the present 
work, each numerical result has been averaged over 10^ 
configurations and convergency has been verified. 


IV. NUMERICAL RESULTS 

In this section we present the non-equilibrium calcu¬ 
lations obtained using the Kwant code^^. An important 
criterion in the investigation of disorder-driven effects is 
the determination of the effective mean free path A of 
the system, as a function of the strength of the disor¬ 
der. To do so, we used the semiclassical formula for the 
conductance 

G = Go/ (l + x) > (9) 

where G is the actual conductance of the sample. Go is 
the ballistic conductance, L is the length of the sample 
and A is the mean free path. This formula indicates that 
as long as the mean free path is well defined, the ratio 
Aeff = I//(l — Go/G) should remain constant when vary¬ 
ing L. In principle, Eq. (9) is only valid as long as lo¬ 
calization effects are weak and when SOC is absent. The 
effective mean free path extracted from our numerical 
calculations are displayed in Fig. 3 (the conductance is 
shown in the inset, for reference) and exhibits no depen¬ 
dence as a function of the length of the sample (except 
for very weak disorder - blue symbols - which indicates 
that the transport is not in the diffusive regime). These 
straight lines constitute a sane signature that the regime 
is diffusive in this parameter space. From now on, we 
choose the disorder strength in the range where the mean 
free path in Fig. 3 shows a constant profile. 



40 60 

Length (a) 


FIG. 3. (Color online) Mean free path as a function of the 
sample length for different disorder strengths. The lower 
curves correspond to the normal (localization free) diffusive 
regime whereas a small divergence is obtained for the weaker 
disorder due to size effects (blue curve). The inset shows the 
corresponding conductances. 


A. Crossover between spin swapping and spin Hall 
effects. 

In order to study the crossover between SHE and SSW, 
we consider the structure depicted in Fig. 2(a). The cur¬ 
rent is injected from the FM layer into the NM layer along 
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FIG. 4. (Color online) Two dimensional mapping of the 
(a,c) X- and (b,d) ^-components of the spin density (in unit 
of 10~^h/2) demonstrating the crossover between SSW and 
SHE, for (a,b) a = 0.3 and (c,d) a = 0.6. The disorder 
strength is T = leV. 


x-direction. Figure 4 displays the two dimensional map 
of the X- and ^-components of the spin density, Sx and 
Sz, for different strengths of the SOC unitless parameter, 
q/= 0.3 and (a=0.6. The spatial profile of the three compo¬ 
nents of the spin density taken along a line passing by the 
center of the sample is shown in Fig. 5. Since the magne¬ 
tization of the FM layer is aligned along the ^-direction, 
the injected current is spin polarized transversally to the 
direction of injection. The transverse component of the 
spin density Sy simply relaxes in the NM due to SOC- 
driven spin relaxation [see Fig. 5, green symbols]. Inter¬ 
estingly, the injected spin current immediately experi¬ 
ences spin swapping, which results in the accumulation of 
the spin density Sx along the edges of the sample, close to 
the FM/NM interface. This accumulation vanishes over 
a length of the order of the spin diffusion length. Simul¬ 
taneously, as the injected spin current gets progressively 
depolarized through spin relaxation, a spin density Sz 
accumulates at the edges of the sample due to SHE. In 
contrast with Sx, Sz survives away from the FM/NM 
and is maintained throughout the length of the sample. 
This crossover between SSW and SHE is therefore con¬ 
trolled by the spin relaxation length, in agreement with 
the scenario exposed in Section II. 

It clearly appears that the injected spin current is 
converted into Jy upon scattering in the normal metal, in 
agreement with the phenomenological prediction of Lif- 
shits and Dyakonov^^ and with the drift-diffusion model 
given by Shchelushkin and Brataas^^. One can recover 
the spin density profile induced by SSW in Figs. 4(a,c) 
by solving the diffusive transport equations, Eqs. (3)-(4), 
to the lowest order in SOC. In the configuration adopted 
in the present work, the spin accumulation injected in 
NM has the form /i^ = where /ig depends 

on the material parameters'^ (conductivity, polarization, 
interfacial resistance etc.). In the presence of SOC, the 
magnitude /io acquires a dependence on y and one finds 
that the spin accumulation at the edges is 

MsU = (2^so/3)/iotanh(L/A5/)e“*/^'’'f. (10) 


This solution shows that the spin relaxation plays an im¬ 
portant role in the build-up of the spin accumulation in¬ 
duced by SSW. In the absence of spin relaxation, no spin 
current is generated and therefore SSW is inactive. But 
on the other hand, when spin relaxation is too strong, ob¬ 
viously, the spin accumulation at the edges of the sample 
does not extend away from the interface. 

(a) q;= 0.3 (b) a=0.6 




FIG. 5. (Golor online) Spatial profile of the {x, y, z)- 
components of the spin density calculated at the upper edge 
of the sample {y = T/2). The blue and red curves represent 
the spatial profile of the spin accumulation induced by SSW 
and SHE, respectively, while the green curve shows for the 
relaxation of the injected spin current. Note that the first 20 
layers compose the ferromagnetic edge. The parameters are 
the same as in Eig. 4. 


B. Dependence on the SOC strength 

The relative magnitude of the spin density accumu¬ 
lated through SSW and SHE depends on the strength 
of the SOC parameter in a nonlinear manner, which is 
at first sight quite surprising since both effects are, in 
principle, proportional to the SOC strength. However, 
the spin relaxation is detrimental to the spin swapping 
since it limits the amount of spin density available for 
the swapping. Indeed, Fig. 6 shows the two dimensional 
mapping of Sx for different strengths of SOC. The in¬ 
duced spin density is localized at the edges of the sample 
and extents over Agf away from the interface. The spin 
relaxation length decreases when the SOC strength a in¬ 
creases from a = 0.2 (Agf ^ 100a) to a = 0.5 (Agf ^ 40a). 

The interplay between the scattering events mentioned 
above in conjunction with the spin relaxation lead to the 
emergence of two regimes which characterize the depen¬ 
dence of the SSW accumulation on the SOC strength a. 
The weak SOC regime corresponds to the enhancement 
of the spin swapping phenomena simply driven by the 
increase of the conversion efficiency with SOC^^, while 
the spin relaxation remains limited. For relatively larger 
values of SOC {a > 0.3), the SSW accumulation de¬ 
creases in magnitude when a increases (Fig. 7). The x- 
component of the spin density varies from ^ 6x 10~^h/2 
to ^ 4 X 10~^h/2 for a varying from 0.3 to 0.5. Indeed, 
while the conversion efficiency increases with a, the spin 
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FIG. 6. (Color online) Two dimensional mapping of the x- 
component of the spin density for various SOC strengths, 
ranging from a = 0.2 to a = 0.5. we choose F = 1 eV in 
all the subplots. 


relaxation reduces the amount of spin current that 
is available for swapping. At large enough a, the spin 
relaxation dominates and the spin swapping decreases. 
When a increases, the spin relaxation length (the extent 
of the spin accumulation at the edges) decreases from 
\sf = for a = 0.2 to Xsf = X^j^/2 for a = 0.4. It 
is remarkable to notice that while SSW has a non linear 
dependence as a function of a, SHE homogeneously in¬ 
creases with a. Indeed spin relaxation, as it depolarizes 
the injected spin current in NM, favors SHE. 



FIG. 7. (Color online) x- and ^-components of the spin den¬ 
sity calculated for one chosen point at the upper edge of the 
sample with the coordinates (xo, T/2). SSW-induced spin ac¬ 
cumulation hrst increases and then decreases when increasing 
the SOC parameter a, while SHE-induced spin accumulation 
monotonously increases. We have verihed that the same be¬ 
havior holds for different values of xq. Here F = 1 eV. 


C. Influence of the disorder 

As mentioned in Section HI, the disorder is a manda¬ 
tory ingredient in our model to get all the effects previ¬ 
ously discussed, a fact that elucidates the pure extrinsic 
origin of the SSW effect. In order to identify the impact 
of disorder on SHE and SSW, we numerically calculate 


the spin density profile for various strengths of disorder 
(i.e. different mean free paths). Similarly to what has 
been noted for the dependence of the SSW accumulation 
on (a, the disorder dependence exhibits the two regimes 
mentioned above. While the SOC strength a controls 
the coupling between the spin and orbital momenta, the 
disorder strength F controls the scattering probability off 
the impurities. 


(a) r=0.6eV (b) SS\N{<5Sx>ilO-^h/2)) 
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FIG. 8. (Color online) Left panels: Two dimensional mapping 
of the x-component of the spin density for (a) F = O.GeV and 

(c) F = leV; Right panels: Spatial prohle of the (b) x- and 

(d) ^-components of the spin density calculated at the upper 
edge of the sample {y = T/2), for different disorder strengths 
(F =0.6, 0.8 and 1 eV). Here, we set a = 0.4. 


Figure 8 represents the spatial profile of the x- and 2 :- 
components of the spin accumulation at the edge of the 
sample for different disorder strengths. In these calcula¬ 
tions, the effective mean free path ranges from A = 100a 
(F = 0.4 eV) to A = 25a (F = 1 eV). Quite remarkably, 
while the extension of the SSW-induced spin accumu¬ 
lation within the normal metal decreases when increas¬ 
ing the disorder strength [see Figs. 8(a) and (c)], the 
overall magnitude of the edge spin accumulation is only 
weakly sensitive to F [see Figs. 8(b)]. This indicates that 
the magnitude of the SSW effect is mostly controlled by 
the spin relaxation strength, while the SSW efficiency 
itself remains independent on the disorder. This obser¬ 
vation is consistent with the diffusive model which states 
that the SSW-induced spin current is proportional to the 
SOC strength only [see last term in Eq. (4)], while the 
edge spin accumulation is controlled by the spin diffusion 
length [see Eq. (10)]. This behavior sharply contrasts 
with the dependence on the SHE-induced edge spin ac¬ 
cumulation displayed in Eig. 8(d). Indeed, our calcula¬ 
tions show that the magnitude of the SHE-induced spin 
accumulation is on contrary very sensitive to the disor¬ 
der strength and almost vanishes for F < 0.4eV, whereas 
SSW survives. This is, again, consistent with the diffu¬ 
sive model, Eq. (4), that shows that in our model ex¬ 
trinsic SHE necessitates a large amount of disorder to 
emerge. 
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V. OBSERVING SPIN SWAPPING 


Before concluding this study, we propose three setups 
for the experimental detection of the spin swapping. The 
first setup, displayed on Fig. 9(a), relies on direct spin in¬ 
jection. A spin-polarized charge current is injected along 
the direction z from a polarizer (denoted FI) into a nor¬ 
mal metal possessing strong SOC (red layer). Since the 
magnetization of the polarizer FI is aligned along 2 ;, the 
spin swapping generates a spin current flowing along the 
direction x and spin-polarized along x. This spin current 
accumulates spins (red arrows) at the interface with a 
second ferromagnet, F2. By sweeping the magnetization 
of F2 along (or by simply reversing the direction of 
the injected current), one should collect a magnetoresis¬ 
tive signal arising from the interfacial spin accumulation 
produced by the spin swapping. The main constraint of 
this setup is that the area of the ferromagnet F2 in con¬ 
tact with the normal metal must be at most of the order 
of Ag£ to maximize the detection. 

The second method we propose, displayed on Fig. 
9(b), relies on non-local spin injection and may be sim¬ 
pler to fabricate. A spin-polarized charge current is in¬ 
jected along —cc, from a polarizer (denoted FI) into the 
normal metal (red layer). A pure spin current, free of 
charge and spin-polarized along the magnetization di¬ 
rection of FI aligned along 2 :, diffuses along (blue 
arrows). Through spin swapping, a spin current polar¬ 
ized along X and flowing along 2 : is injected into a sec¬ 
ond ferromagnet F2 placed on top of the normal metal. 
Again, by sweeping the magnetization of F2 along 
or by reversing the direction of the injected current, one 
should collect a magnetoresistive signal arising from the 
spin swapping. The main constraint of this setup is that 
the distance between FI and F2 must be smaller than 
Agf. This setup presents some similarities with the fa¬ 
mous spin field-effect transistor proposed by Datta and 
Das^^, whose operation relies on coherent spin precession 
induced by Rashba SOC rather than spin-orbit coupled 
impurities. 

The two previous setups present the inconvenience of 
being limited by the spin diffusion length of the spin-orbit 
coupled metal. One alternative situation is depicted in 
Fig. 9(c), which is a variant of the setup displayed in Fig. 
9(a). This structure consists in a Hall cross where a spin 
current is injected from a first ferromagnet FI into the 
cross and spin-orbit coupled impurities are only present 
in the center of the cross. For instance, if one imagines 
a Hall cross made of Cu, doping the center of the cross 
with, say. Bismuth impurities would allow for the spin 
swapping to take place at this location only while pre¬ 
venting spin relaxation in the branches of the cross. The 
detection of the spin swapping-induced spin current can 
be achieved by performing a magnetoresistance measure¬ 
ment on F2. Another advantage is that the magnetiza¬ 
tion direction of FI and F2 are both in-plane. 



FIG. 9. (Color online) Sketches of various setups for mea¬ 
suring spin swapping: (a) direct injection using a three ter¬ 
minal devices, (b) non-local injection in a one-dimensional 
’’spin transistor” configuration and (c) direct injection in a 
Hall cross. The blue layers represent ferromagnets and the 
red layer represents a normal metal. The thick red arrow 
denotes the direction of injection while the small blue (red) 
arrows indicate the direction of the electron spins before (af¬ 
ter) spin swapping. 


VI. CONCLUSION 


We have shown that spin-orbit coupled disorder in¬ 
duces a strong SSW that may even exceed extrinsic SHE 
in metals. The crossover between these two effects is 
controlled by the spin relaxation of the injected spin cur¬ 
rent. While both SSW and SHE emerge through scat¬ 
tering on spin-orbit coupled disorder, these two effects 
present a noticeably different dependence as a function of 
both the SOC and disorder strengths. Indeed, while SHE 
monotonously increases with the SOC strength, SSW dis¬ 
plays a more complex behavior due to the competition 
with spin relaxation. In addition, while SHE is very sen¬ 
sitive to the disorder strength, it turns out that the SSW 
efficiency is rather controlled by spin relaxation than by 
the disorder strength itself. 

Einally, we propose three experimental setups that 
shall allow for the experimental detection of the spin 
swapping in metals possessing impurities-driven spin- 
orbit coupling. Materials in which extrinsic SHE dom¬ 
inates over the intrinsic contributions, such as heavy 
metal-doped light metals such Au(Pt)^^, Au(W)^^, 
Cu(Ir)^^, Cu(Bi)^^, or even rare-earth doped metals^^, 
are attractive candidates for the observation of this effect 
as they associate strong spin Hall angle with a reasonably 
long spin diffusion length. 
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